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O ■ Abstract 

In principle, the complete spectrum and bound-state wave functions of a quantum 
field theory can be determined by finding the eigenvalues and eigensolutions of its 
light-cone Hamiltonian. One of the challenges in obtaining nonperturbative solu- 
tions for gauge theories such as QCD using light-cone Hamiltonian methods is to 
renormalize the theory while preserving Lorentz symmetries and gauge invariance. 
For example, the truncation of the light-cone Fock space leads to uncompensated 
ultraviolet divergences. We present two methods for consistently regularizing light- 
cone-quantized gauge theories in Feynman and light-cone gauges: (1) the intro- 
duction of a spectrum of Pauli-Villars fields which produces a finite theory while 
preserving Lorentz invariance; (2) the augmentation of the gauge-theory Lagrangian 
with higher derivatives. In the latter case, which is applicable to light-cone gauge 
(A + = 0), the A~ component of the gauge field is maintained as an independent 
degree of freedom rather than a constraint. Finite- mass Pauli-Villars regulators can 
also be used to compensate for neglected higher Fock states. As a test case, we apply 
these regularization procedures to an approximate nonperturbative computation of 
the anomalous magnetic moment of the electron in QED as a first attempt to meet 
Feynman 's famous challenge. 
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1 Introduction 



The gyromagnetic ratio of the electron g e - = 2.0023193043768(86), the ratio 
of the spin precession frequency to the Larmor precession frequency in a static 
magnetic field, is an intrinsic property of an individual lepton. It is now known 
experimentally to 12 significant figures [1] - the most precisely known funda- 
mental physical parameter. The anomalous moment a e = ^—, the deviation 
of the gyromagnetic ratio from Dirac's value g e = 2 due to quantum fluc- 
tuations, has now been evaluated through order a 4 in perturbative quantum 
electrodynamics [2,3]. 

At the 12th Solvay Conference, Feynman presented a challenge [4]: "Is there 
any method of computing the anomalous moment of the electron which, on 
first approximation, gives a fair approximation to the a term and a crude one 
to a 2 ; and when improved, increases the accuracy of the a 2 term, yielding a 
rough estimate to a 3 and beyond." An interesting attempt to answer Feyn- 
man's challenge using sidewise dispersion relations was pioneered by Drell and 
Pagels [5], but it is difficult to make this method systematic. 

The anomalous moment of a spin-half particle can be evaluated without ap- 
proximation from the overlap of its light-cone Fock-state wave functions. The 
overlap of the two-particle one-fermion-one-boson light-cone Fock state I/7 > 
yields Schwinger's contribution [6] a e = f-. The light-cone wave functions 
with n > 2 QED quanta contribute to a e beginning at order a n ~ x as well as 
higher orders. Thus a systematic evaluation of the lepton's Fock-state wave 
functions as an expansion in Fock number rather than perturbation theory 
would provide a physically appealing answer to Feynman's challenge [7]. 

In principle the complete spectrum of a quantum field theory can be deter- 
mined by finding the eigenvalues of the light-cone Hamiltonian [8]. The Fock- 
state expansion of the eigensolutions at fixed light-cone time r = x + = x° + x z 
provides a frame-independent wave-function description of the elementary and 
composite states in terms of the quanta of the free Hamiltonian. The dis- 
cretized light-cone quantization (DLCQ) method utilizes periodic boundary 
conditions to truncate the size of the Fock-state expansion while preserving 
boost invariance. This method has been successfully applied to a large variety 
of gauge theories and supersymmetric theories in 1 + 1 and 2 + 1 dimensions. 

An essential problem in applying light-cone Fock-state methods to renormaliz- 
able gauge theory is to regularize the calculations in such a way as to preserve 
Lorentz and gauge invariance, or at least preserve them well enough to allow 
an effective renormalization to be performed. Since the method of regulariza- 
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tion must also allow for efficient calculations to be performed, the problem 
presents a challenge. 

We have recently performed nonperturbative calculations using the generalized 
Pauli-Villars (PV) method as an ultraviolet regulator of (3+l)-dimensional 
quantum field theories. We include a sufficient number of PV fields in the 
Lagrangian to ensure that perturbation theory is finite. This method explic- 
itly preserves Lorentz invariance; in some cases, such as QED, it effectively 
preserves gauge invariance. In a case where it breaks gauge invariance, such 
as QCD, we have to add counterterms. The PV regularization can produce a 
finite theory which preserves Lorentz and gauge symmetries. However, if we 
do not have the exact solution, we must develop approximate methods. Our 
approximation involves truncating the Fock space. The truncation will break 
all of the symmetries. However, the usefulness of the truncated answer is a 
question of accuracy rather than an issue of symmetry breaking. With reg- 
ulators in place we presume that there is an exact solution which preserves 
all symmetries including gauge invariance, and if our approximate solution is 
close to the exact one, even if the small difference is in such a direction as to 
maximally violate the symmetries, it is still a small difference. Of course, the 
inclusion of negative-metric fields in the Lagrangian will also violate unitarity. 
We shall have more to say about this issue below. 

The generalized PV method has been applied successfully to Yukawa-like the- 
ories [9,10,11,12,13] where there are no infrared divergences and no need to 
protect gauge symmetry. An important conclusion of these studies is that past 
some threshold (which depends on the values of the coupling constant and the 
values of the PV masses), there is always a rapid drop off of the projection of 
the eigensolution wave function onto higher Fock sectors, in contrast to the 
equal-time Fock-space expansion. This provides a strong motivation for the 
light-cone representation as a viable approximant to nonperturbative theory. 

In this paper we will test the convergence of the PV-regulated light-cone Fock- 
state expansion for (3+l)-dimensional gauge theory by applying it to a non- 
perturbative calculation of the electron anomalous moment in QED. In some 
ways this application to QED is not an ideal test of these nonperturbative 
methods: the physical electron is a very perturbative object. We do not ex- 
pect to do better, or even as well as perturbation theory. However, that is not 
our objective; we simply want to verify that an approximate nonperturbative 
solution for the electron's magnetic moment is an approximation to QED. 
Somewhat related work, from a strictly perturbative point of view, can be 
found in [14]. 

Careful studies have also shown [15,16] that the perturbative series obtained 
from particular combinations of PV fields and higher derivative regulators 
give the same result as standard perturbation series regulated with dimen- 



sional regularization. These studies not only included Yukawa theory, but also 
non-Abelian gauge theory. In Section 4, we shall apply this method of reg- 
ularization, together with a Fock-space truncation to the calculation of the 
electron magnetic moment. This is the first application of this regularization 
to a nonperturbative problem. 

We find that there are three problems which must be solved in order to pro- 
duce a useful calculation of the electron's magnetic moment: the problem of 
uncanceled divergences, the problem of maintaining gauge invariance, and the 
problem of new singularities. We believe that we have found effective solutions 
to these problems, at least for the present calculations. The problem of un- 
canceled divergences occurs anytime we truncate the Fock space. For example, 
if we truncate the physical electron's Fock space to include only the subspace 
of one fermion and one photon, calculate the wave function nonperturbatively, 
and use that wave function to calculate the moment, we obtain a result of the 
form: 



a 



x [finite quantity] 



1 + a x [finite quantity] + a x [finite quantity] log 



Til' 



where fi± is the PV photon mass. If we let /xi become infinite, we will obtain 
a zero anomalous moment. 

The origin of the uncanceled divergence in Eq. (1) can be seen by examining 
the two-loop contributions to the electron moment in perturbation theory. 
The relevant double-ladder Feynman diagrams are shown in Fig. 1. The Dirac 
Fi(q 2 ) and Pauli i 7 2(? 2 ) form factors correspond to S z = ±S' Z matrix elements, 
respectively. (The q + = 0, q\ = Q 2 = —q 2 frame is assumed.) Fock states with 
particle number 1,2,3 contribute as indicated in the figure. Only time orderings 
with positive k + appear in light-cone quantization. The Ward identity Z\ = Z 2 
guarantees the cancellation of the divergences from wave function and vertex 
renormalization subgraphs. However, if the three-particle Fock state I/77 > is 
excluded by the truncation of the Fock space, F 2 is ultraviolet (UV) divergent 
since the vertex correction shown in amplitude (c) is missing. The divergent 
contribution to the lepton anomalous moment at = ^(O) is the source of the 
UV divergence in the denominator of Eq. (1) due to Fock-space truncation. 
Of course, iq remains UV finite. The divergence in Eq. (1) does not happen 
in perturbation theory; since the numerator is already of order a, we would 
use only the 1 from the denominator. At order a 2 there would be new terms 
in the numerator which would cancel the divergent terms in the denominator 
of Eq. (1). 

The above discussion illustrates the problem of uncanceled divergences. While 
we could find ways of allowing the bare mass and the coupling constant to 
depend on the PV masses to give a finite expression, the results would not look 
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Fig. 1. Light-cone time-ordered contributions to the lepton form factors, corre- 
sponding to the order-a 2 ladder Feynman diagram in perturbative QED. The verti- 
cal dashed lines mark intermediate states with the indicated number of constituents. 
For amplitudes associated with diagrams in (d), with loops on the external legs, a 
factor of 1/2 is applied. 



anything like the results from perturbation theory, since in QED the coupling 
is renormalized only by vacuum polarization. In addition, the results would 
not make sense physically. Our resolution of this difficulty is to keep the PV 
masses finite. The motivation is as follows: If the limit of infinite PV masses 
would give a useful answer in the case where we do not truncate the Fock 
space (so we have no uncanceled divergences), then there must be some finite 
value of the PV masses that would also give a useful answer. The question is 
whether we can use a sufficiently large value. To answer that question we must 
consider that there are two types of error associated with the values of the 
PV masses. The first type of error results in having these masses too small; 
then our wave function will contain too much of the negative-normed states, 
unitarity will be badly violated, and in the worst case we might get negative 
probabilities. We can roughly estimate the magnitude of that type of error as 

E 1 = M 2 jMl (2) 



where Mq is the physical mass scale and M\ is the PV mass scale. The other 
type of error results when the PV masses are too large; in that case the true 
wave function will project significantly onto the parts of the representation 
space excluded by the truncation. We can roughly estimate the magnitude of 
that type of error as 

E2 = m^ry (3) 



where |<&+) is the projection of the wave function onto the excluded sectors. In 
practice, the projection onto the first excluded Fock sector can be estimated 
perturbatively using the projection of P~ onto the higher sectors as the per- 
turbing operator. Without additional information, the best that we can do is 
to set Ei equal to the perturbative estimate of E 2 . Below, we will apply this 
procedure to the calculation of the magnetic moment done in this paper. In 
general, if, at our estimated optimum value for the PV mass scale both types 
of error are small, we can do a useful calculation; otherwise, we cannot do a 
useful calculation without expanding the part of the representation space that 
we include in our calculation. 

The main reason for believing that we can do a useful calculation in spite of 
the problem of uncanceled divergences is the lesson from the earlier studies 
mentioned above: the observed rapid drop off of the projection of the wave 
function onto higher Fock sectors. Just where this rapid drop off occurs de- 
pends on the theory, the coupling constant, and the values of the PV masses. 
At weak coupling and relatively light PV masses, only the lowest Fock sectors 
are significantly populated. At stronger coupling or heavier PV masses, more 
Fock sectors will be populated; but eventually the projection onto higher sec- 



tors will fall rapidly. The rapid drop off of the projection of the wave function 
onto sufficiently high Fock sectors is the most important reason why we do 
our calculations in the light-cone representation. For any practical calcula- 
tion in a realistic theory, we have to truncate the space, and we must have 
a framework in which that procedure can lead to a useful calculation. The 
rapid drop off in the projection of the wave function will not happen in the 
equal-time representation, mostly due to the complexity of the vacuum in that 
representation. 

These features can be explicitly demonstrated by setting the PV masses equal 
to the physical masses. In that case the theory becomes exactly solvable [12]. 
The spectrum is the free spectrum, and the theory is not useful for describ- 
ing real physical processes due to the strong presence of the negative-norm 
states in physical wave functions; however, it still illustrates the points we 
have been trying to make. In the equal-mass case, the physical vacuum is the 
bare light-cone vacuum, while it is a very complicated state in the equal-time 
representation; the physical wave functions project onto a finite number of 
Fock sectors in the light-cone representation but onto an infinite number of 
sectors in the equal-time representation. 

As the PV masses become larger than the physical masses, the light-cone wave 
functions project onto more of the representation space. This effect is increased 
as the coupling constant becomes larger. However, the wave functions remain 
much simpler than in the equal-time representation, and, to the extent we can 
do the calculations, there is always a point of rapid drop off of the projection 
onto higher Fock sectors. Due to this rapid drop off, we expect to find a PV 
mass scale such that the error in the calculated value for a given observable 
from the presence of negative-norm states and the error from truncation can 
both be made arbitrarily small. Therefore, we believe the requirement to keep 
the value of the P V masses finite does not impose a limit on the accuracy which 
could, in principle, be achieved. In practice, the size of the representation space 
may be too large for presently available computing facilities. Furthermore, for 
the method to be useful in practice, there must not only be a value of the PV 
mass for which both types of errors are small, but there must be a wide range 
of such values since the optimum value for the PV mass can only be rather 
crudely estimated. A principal objective of the present work is to test these 
ideas on a physically realistic problem to which we know the answer. 

The problem of maintaining gauge invariance turns out to be nontrivial and 
quite instructive. In the next section we consider the following procedure: 
we take the light-cone evolution operator P~ for QED in light-cone gauge 
A + = as constructed, for example, in Refs. [17] and regulate it with PV 
fields. We looked at several cases: PV photons alone, PV fermions alone, and a 
combination of both. This procedure does not lead to a viable nonperturbative 
formalism for the electron magnetic moment; furthermore, it does not lead to 



a correct calculation of the electron self-energy, even at order a. The problem 
can be traced to a failure to maintain gauge invariance. The breakdown of the 
straightforward implementation of the PV method in light-cone gauge shows 
that the successful construction of the nonperturbative theory is nontrivial. 

In Section 3 we perform the calculation in Feynman gauge with one PV photon 
and one PV fermion; this leads to a consistent formulation for the nonpertur- 
bative calculation of the electron moment in QED. In Section 4 we perform 
the calculation in light-cone gauge, but use the more sophisticated method 
of regulation proposed in [16]. We show that this method provides a success- 
ful formulation of the nonperturbative electron moment problem in light-cone 
gauge with a result very similar to that in Feynman gauge. 

In Sections 3 and 4 we face the problem of new singularities. We have to do 
integrals with denominators of the form (— M 2 x(l — x) +m 2 x + fi 2 (1 — x) + k 2 L ), 
where M is the physical electron mass, m is the bare electron mass, and [i 
is the photon mass. When the bare mass is less than the physical mass, as 
is the case in QED, there can be a zero in this denominator. In perturbation 
theory the expansion is about M = m, and the denominator cannot vanish as 
long as the photon is given a small nonzero mass. The standard techniques in 
perturbation theory thus avoid this singularity. We find that when the zero is 
a simple pole, the principal-value prescription is correct. However, in the wave 
function normalization the denominator is squared, so there is a double pole, 
and we must give it a meaning. We propose the following prescription: 



jdy 



dk\ 



f(y, k] 



m 2 y + f i 2 (l-y)-M 2 y(l-y) + k 2 ± } 2 
J|m — J dy J dk 2 ± f(y,k 



[m 2 y + ^§(1 -y)- M 2 y{\ - y) + k\ - e] 

(4) 



[m 2 y + fi 2 (l -y)- M 2 y{\ -y) + k\ + e 

where simple poles are prescribed as principal values. 

This prescription has the interesting consequence that the wave function nor- 
malization is infrared finite whereas it is infrared divergent in perturbation 
theory. Given this prescription, the true singularity occurs at M = m + /i; 
in perturbation theory, with M = m, this is at \x — 0, which is the infrared 
singularity and the reason that the photon mass cannot be taken all the way 
to zero in perturbation theory. For the nonperturbative calculation, the phys- 
ical photon mass can be taken to zero since M ^ m. The basic requirement 
of these prescriptions is that they preserve the Ward identities. We antic- 
ipate using this prescription for QCD where the basic requirement will be 
the preservation of the Ward-Takahashi identities. We have not shown that 



the prescription preserves the Ward identities in QED, but it does lead to a 
successful calculation in the present case. 

A different approach to this same dressed-electron problem has been taken by 
Karmanov, Mathiot, and Smirnov [18]. They use a covariant form of light- 
cone quantization without Pauli-Villars regularization. Their Hamiltonian 
then contains instantaneous fermion interactions, and, in Feynman gauge, the 
infinite number of terms generated by inversion of the covariant derivative. 
The problem of uncancelled divergences is avoided by application of sector- 
dependent renormalization [19]. They must construct counterterms explicitly. 
However, they truncate in a Fock basis where the constituent electron has 
the same mass as the dressed electron, bringing their calculation closer to 
perturbation theory; in fact, they find that any signifant difference with per- 
turbation theory will not appear until the basis is expanded to include higher 
Fock states. Also, they do not calculate the anomalous moment. 



2 Trouble in Light-Cone Gauge 

In this and the next section the notation that we use for light-cone coordinates 

is 

x ± = x° ± x 3 , xj_ = (x 1 , x 2 ). (5) 

The time coordinate is x + , and the dot product of two four- vectors is 

p • x — -(p + x~ + p~x + ) — pj_ ■ x±. (6) 

The momentum component conjugate to x~ is p + , and the light-cone energy 
is p~ . Light-cone three- vectors are identified by underscores, such as 

p = (p + ,p±). (7) 

We use the following choice for the 7 matrices 





For additional details, see Appendix A of Ref. [9]. 

We use the standard P~ for light-cone quantized QED in light-cone gauge [17] 
with modifications due to the inclusion of the PV fields. We should remark, 
however, that with the inclusion of any number of PV Fermi fields, the four- 
point interactions which would take a state of one electron and one photon to 
another state of one electron and one photon are missing from P~; that such 
terms are not included below is not an omission; the calculation is complete 
in our chosen subspace. We truncate the Fock space to the one-fermion sector 
plus the one-fermion, one-photon sector. We then solve the eigenvalue problem 

P+p-\s) = M 2 |s), (9) 



where the total P± of the state is taken equal to 0. As always with a Tamm- 
Dancoff truncation, we can solve for the wave function in the highest Fock 
sector (one fermion plus one photon) by hand and obtain an equation in the 
one-fermion sector. 

We have regulated the theory in several ways. We shall describe one particular 
choice in some detail: We use three PV Fermi fields with flavor-changing cur- 
rents. The flavor-changing currents break gauge invariance, and thus we might 
expect to require counterterms to correct for that. Nevertheless, we will pro- 
ceed with the calculation without counterterms. That will allow us to study a 
case where a proper respect for gauge invariance is not maintained. Also, as 
we shall argue below, the source of the breaking of gauge invariance is much 
deeper than that due to the flavor- changing currents. 

For the Lagrangian we take 

-F^F, V + ]T -&(ird» - rmW - e^W4*> ( 10 ) 



4 p " ^ v 

^ i=0 "* 



where 



3 3 3 3 

i=0 i=0 i=0 i=0 



A particular realization of these PV conditions is 

{vi} = {1; 3; -1; -3}, {m^} = {m ; -m + 2m 3 ; -2m + 3m 3 ; m 3 }.(12) 

With this choice it is convenient to define 

V'o = fp'o, ipi = VSip[, ip 2 = ip' 2 , ip 3 = VSip' 3 , (13) 
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so that the ip fields are canonically normalized (except for the minus signs for 
ip2 and ■0 3 ). The coupling to the A field is 



e^WA* = J29iA^3 A ^ 



(14) 



where 



#00 = #02 = #22 = e, 5-01 = #03 = #12 = #23 = v^e, 
#11 = #13 = #33 = 3e, 

and ftj = ^j. 

We use the mode expansions 
1 



(15) 



A+(x) 



A\x) 



V167T 3 



Y, d kXs b s (i, k)e-** + dl s (i, k)e +i ^ 



VT67T 3 x 



]T / dk -j= a A (fc)e i (A)e-^ + aUi)e i (A) e 



-ik-x 



(16) 
(17) 



where the polarization states are 



X+i 



V2 





1 



X-h 



V2 



1 


v-iy 



e±,+i = 77|(M) 



1 



^,-1 = ^11,-) 



(19) 



The vertex functions can be computed from P as 
Vij = VlL^ VI ~ x x(l - x) {kx ~ iKl 



^v ~ 7t^^~r == ^.T\ _ x y ^ x ^y' : 



V167T 3 \J\ — X 



u-- 6ij ' 



m~ 



I — TYl ■ -I- 

VT6^VT=x K * x 



77.. - _!i! 1 



W; 



, ( —mi + 



(20) 
(21) 
(22) 
(23) 
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where 



e .. = j9ij, if i = 0, 1; 

e ^-i-^, if i = 2, 3. [ZV 



We expand the eigenstate as 

3 3 . 

\s >=J2zib\.(i,P)\0 > + E / dkf(i,x,k ± )bl(i,k)al(P-k)\0 > 

j=0 i=0 

3 — /• 

+ E E dkg s (t,x,k ± )bl(t,k)al s (P-k)\0> . (25) 

We will take the wave function normalization zq to be 1 for the moment and 
calculate it later. From P + P~\s >= m 2 e \s > we find immediately (in units 
where m e is taken to 1): 

Hi r k,)- ^ 3 3=o z 3 V 3i( x > k -J (2 a\ 

x l—x 



E?=o ZjWjfak 

m 2 +k 2 fj, 2 +k 2 ' 



g+(i, x, k±) = J \]J l \, h 2 , (27) 



l-x 



Z^i=o ZjU jiix, k_i 

rn^+k 2 ^ ^ 2 +k 2 



9-&*m) = 3 z,lj ;, a r; a ■ (28) 



l-x 



From all this we derive the four nonlinear equations 

V* i (x,k ± )Ef=oZiV lj (x,k±) 



Zim%+ ^2 / dxd 2 kj 



3=0 1 „ 1 



x l—x 



+ > / arret A; j. — - 2TT2 n — t 

s l—x 



: ' '■ Uji(x,k ± ) Ef= *j #y (&, *. 

i=o •> 1 - 



E/^x^^f^P = ,, (29) 



z l—x 

These can be written in terms of the three integrals: 



J = dxdz^ — r-i — '-j- r , (30) 

J ~r, x x(l — x — mf(l — x — a z x — z 



i=0 
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I 1= dxdz^-r- r 2 * - , (31) 

J ~q x{± — x) — mf{l — x) — fj, z x — z 



'' XV; 



h = / dxdz J2 —pi \ ITT S 2 ' ( 32 ) 

J r, x(\ — x) — rriid — x) — ux — z 



., v j x) — raf (1 — x) — n 2 x — z ' 
where z = k\. The nonlinear equations become 



e 2 3 e 2 3 

iU7r j=o iU7r j=0 

e 2 3 
+ — —I V ejiZjinifrij = z { . (33) 

167T 2 ~^ 

While these equations are complicated, they can be simplified by the observa- 
tion that, for the parameter values of interest, J is very much larger than I 
and I\. Also, for those parameter values, zi, Z2, and z% are small compared to 
one, and a solution sufficiently accurate for our needs is given by the simple 
expression 

ml = 1 - ^ J- (34) 



We can use the wave function to calculate the anomalous magnetic moment 
using the formalism of Brodsky and Drell [20]. The contribution from the 
physical field is much larger that the contributions from the PV fields, so we 
have 

a 2 f i 17liX 2 (l — x) 

a e = ~ z o / dx 2 ; T7~i \ m o (35) 

n J m(x + jJL Z (l — x) — x(l — X) 



where zq is determined by wave function normalization: 

l/z 2 = l + z 2 -zt-z 2 3 (36) 

+ |rfxrf 2 A; ± (|/(0,x,^)| 2 + |/(l,x,^)| 2 -|/(2,x,^)| 2 -|/(3,a:,^)| 2 ) 

+ (\g40,x,k ± )\ 2 + \g + (l,x,k ± )\ 2 -\g + (2,x,k ± )\ 2 -\g + (3,x,k ± )\ 2 ) 
+ (|(?-(0,a;,^)| 2 + |(7_(l,a:,^)| 2 -|( ? _(2,x,^)| 2 -|( ? „(3,x,^)| 2 ). 



We find 
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We can now state the problem with this calculation: J has a very strong 
dependence on m 3 (the PV mass scale). That, in turn, gives our estimate of 
the anomalous magnetic moment a very strong dependence on m 3 . If we use 
units of 7T, so that the correct value is near one, then we find that, even with 
a value for the photon mass as large as 0.5 electron masses, when m 3 changes 
from 3 times the electron mass to 7 times the electron mass, a e changes from 
1.2 to -1.2. If we use a smaller value for the photon mass, which we would 
surely have to do to get useful results, the dependence is even stronger. Since 
we cannot hope to estimate the optimum value for the PV mass scale even 
to within this range, the present calculation is clearly useless. The problem 
is clearly the loss of gauge invariance; gauge invariance should prevent such 
strong behavior. One might reasonably think that the problem is the flavor- 
changing currents we have included in the calculation, but we shall now argue 
that the worst breaking of gauge invariance has a more fundamental source. 
We shall return to the breaking due to flavor-changing currents in the next 
section. 

We note that if we keep only the physical field and set M = mo = m, the 
function which appears in our nonlinear equations is just the (unregulated) 
one-loop fermion self-energy 

a a r 1 tt^g z + m 2 (l - x) 2 

f (J - 2 h + I ) = f dxdz'- - ^ / / 2 .(38) 

Air An J x m z x{± — x) — m z [\ — x) — n z x — z 

Therefore, a very useful point of comparison is the evaluation of the fermion 
self-energy calculated in the paper by Brodsky, Roskies and Suaya [21], here- 
after referenced as BRS. They evaluated all the graphs needed to calculate 
the electron's magnetic moment in perturbation theory through order a 2 . In- 
cluded in their calculations is the one-loop electron self-energy They did not 
use light-cone quantization but wrote down time-ordered perturbation theory 
in the equal-time representation and then boosted to the infinite momentum 
frame. They worked in Feynman gauge, but the electron self-energy should be 
gauge invariant. They give the self-energy as the sum of two separate pieces, 
m a and m&; m& results from boosting the z-graph while m a comes from the 
other time ordering. They obtain in our notation (see Eqs. (3.40) and (3.41) 
of BRS) 

2,1 „ , ™2/i _ Arr . _ rY .2\ 



e 2 r 1 z + m 2 {\ — Ax 



x 



ma = 7^ / dxdz TT\ S 271 S 2 ' ( 39 ) 

lbir z J x m z x{± — x) — m z {± — x) — \i z x — z 



e 2 



m h 



16tt 2 



f dzln(fi 2 + z). (40) 



BRS show that if the theory is regulated with the inclusion of two PV photons, 
the sum of m a and rrib is equal to the usual Feynman one-loop self-energy in 
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QED, regulated in the same way. Notice that ra& implements the requirement 
of chiral symmetry: the shift in the bare mass is zero if the bare mass is zero, 
that is, Srrib = —5m a | m=0 . Terms such as m h are often missed in light-cone 
quantization. That fact has been noticed at least as far back as the work of 
Chang and Yan [22]. Those authors suggest that the problem can be solved 
by including an extra PV field and using it to implement the chiral symmetry 
condition; that is a technique we have used in the past [9] . Since m& does not 
depend on the mass, inclusion of PV fermi fields will also solve the problem. 
That possibility was noticed by BRS, and we shall use this method in the next 
sections. Therefore without PV fermi fields we might not expect to get the 
sum of m a and m& correctly, but we should expect to at least reproduce m a . 
We do not have to get the same integrand, but we should get the same result 
after integration if the regulation preserves covariance. However, if we use two 
PV photons, or three PV photons with the third field used to eliminate mj,, 
we do not obtain the correct result for m a . 

We can gain some further insight into what is going on by observing that, in 
Feynman methods, the light-cone gauge is obtained from Feynman gauge by 
the replacement 

n ■ k 



From that replacement there is no obvious source for the double pole at x — 1 
that we see in our light-cone calculation. The explanation is that, in some 
formal sense without worrying about regulation, light-cone quantization ac- 
complishes the replacement by writing 

( uv nv _ n ^ + nvK , fcVv; _ fcW 

g g n-k (n-k) 2 (n-k) 2 ' l } 



The first three terms correspond to the usual three-point interactions that one 
obtains by usual light-cone quantization as given, for instance, in [17]. The last 
term is given by the so called "instantaneous" (four-point) interactions which 
result from solving the constraint equation for the photon. In at least some 
cases, for tree level processes the required cancellations actually occur, and 
results equivalent to the equal-time formulation are obtained. However, here, 
at one loop, the cancellations are not working correctly, and we can see why 
as follows: If we write down the Feynman integral whose numerator is 

k 2 n^n v , , 

- j^w (43) 
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regulate the calculation with two PV photons, and perform the k integral, 
we get an amplitude which cannot be obtained from a four-point interaction 



i 4-^ 



- m 2 x 



dxdz- 2 —; (1 -*> a ' - 2 . (44) 

x m z x[\ — x) — m z (l — x) — jj, z x 



When this amplitude is added to our light-cone amplitude, (38), we get pre- 
cisely m a . There is no doubt that m a -- with the perturbative denominator, 
m 2 x(l — x) — m 2 (l — x) — fi 2 x — z, replaced by the nonperturbative denomi- 
nator, M 2 x(l — x) — m 2 (l — x) — fi 2 x — z — is the function that should enter 
our nonlinear equations. We could fix the problem here by just using m a with 
the perturbative denominator replaced by the nonperturbative denominator, 
but since we do not know how to make similar corrections in other cases, we 
do not consider that replacement to be useful. It is clear that the problem is 
that gauge invariance has been lost in solving the constraint equation 

d 2 _ A~ + d.diA* = -e^l^+. (45) 



It is possible that the wrong boundary conditions have been used in solving 
this equation or that the equation must be modified: the constraint equation 
satisfied by the regulated fields and something like Schwinger terms may need 
to be included. The loss of gauge invariance using standard light-cone tech- 
niques deserves further study. We will now turn our attention to the use of 
other gauges and other methods of regulation. 

In the next section we will discuss light-cone-gauge quantization using Feyn- 
man gauge and PV regularization. In Sec. 4 we shall show that a successful 
calculation can also be made in light-cone gauge, if the formalism is aug- 
mented with higher derivative regulators and several Pauli-Villars fields. In 
that case, due to the higher derivatives, A" is a degree of freedom, and there 
is no equivalent of Eq. (45). 



3 Feynman Gauge 



In this section we shall calculate the electron's magnetic moment using light- 
cone quantization in Feynman gauge. We shall regulate the theory by the use 
of one PV photon and one PV fermion with the inclusion of flavor-changing 
currents. The Lagrangian is thus 



J2 (-h-lYFFFi^ + (-1)^(17"^ - rmWi + B t d,A? + l -B t B t 
i=o v 4 z 
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-e^-fM" ( 46 ) 

where 

A" = J2A*, V = I>, Fr = PA?-d v AV. (47) 

i=0 i=0 



Here, % — indicates the physical fields, and i = 1, the PV (negative- metric) 
fields. 

We will now discuss two important consequences of including PV Fermi fields 
with flavor-changing currents, one good effect and one apparently bad effect. 
The good effect pertains to the operator P~. If one works out P~ including 
only the physical fields, one encounters the need to invert the covariant deriva- 
tive [23] d- — eA_. The same problem occurs in any gauge where A- is not 
zero. This complication is perhaps the main reason that gauges other than 
light-cone gauge have received relatively little attention in the light-cone rep- 
resentation. While the inverse of the covariant derivative can be approximately 
defined by a power series in e, or, in a truncated space may be calculated ex- 
actly if the truncation is sufficiently severe, it is not clear that P~ has been 
fully specified. However, with the inclusion of the PV fermions with flavor- 
changing currents, this problem does not occur: the inverse of the covariant 
derivative is replaced by the inverse of the ordinary derivative. The part of 
P~ that we shall need in our calculations is given by 









p 
" l + k k-^af(k)at(k) 



+ E / d R d l{ b i,s(p) [ b jAl) V ij,2s(P,g) 
i,j,l,s,fi 

+ bj^siqjUtj^siBH)} 4u(i ~p) + h.c.} 
where e^ = (—1, 1, 1, 1) and 



e p±-q±± ip± xq± + m(mj + p + q + 

v ^ i)= vms rtw=F ■ (49) 

V167T 3 p + q + ^/q + -p + 

yi( nn \- e p + (q 1 ±iq 2 ) + q + (p 1 Tip 2 ) , , 

^' ±fe£) "7l6P p+qWd + -p + ' { } 
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J ~- V167T 3 p + q + y/q+ -p+ 



=Fe rrijlp 1 ±ip 2 ) — m i (q 1 ±iq 2 ) 
V167T 3 p + q + Vl + ~P A 



k&s>= ::, '"ZVzr'l — ■ w 



J ~- V167T 3 p+q+ v / q+ -p+ 

,,i , v ±e m i q + -m j p + 

^ fe£) = 7!lP ^vg^F ' (55) 

rr2 / \ ie rn i q + -m j p + 

^(£,£)= Vlfa3p+9+V9+ _ p+ - <»> 



The apparently bad effect of the flavor- changing currents is that they break 
gauge invariance. That would seem to require the inclusion of counterterms in 
the Lagrangian to correct for the symmetry breaking. It turns out that such 
counterterms are not necessary: as we shall see, we can take the limit of the 
PV fermion mass mi — > oo. One might properly worry that there might still 
be residual, finite effects of the necessary counterterms, but the counterterms 
go to zero as powers of m\ while the only divergences we encounter are logs. 
We shall therefore proceed with the calculation using only the Lagrangian 
given above. 

We use mode expansions similar to those of (16) and (17) and expand the 
wave function as 

|V> = ]>>A+Ga|0> + £ JdkC^4s(k)ai(P-k)\0). (57) 



We set the total transverse momentum of the state to zero. We can solve for 

the C's as 



^ + ,i,l\—) ~ „ m 2 + k 2 u 2 + k 2 ' V / 

c , ( m S;(-l)'"*/^a + (fcP) m] 

^-H\Dl — m 2 +k 2 „ 2 +k 2 ' V / 

x l—x 



The eigenvalue equations for Zi become 

(M 2 - m 2 ) Zi =fdx d 2 k ± J2 (-l) i,+J+ ^(P+) V (60) 



X 



V^ + (t,P)V^(k,P) + UZ J+ (k,P)U»* + (k,P) 



M 2 _ m l +k l _ 4 



1-x 



which can be written more usefully as 

(M 2 - m1)zi = 2e 2 ^(-l) j [m 2 Zj J + m^m^o (61) 

3 

— 2Mzj(mi + mj)Ii , 
with 



dxdkj v (-1)'+' {m 3 /M) n 

16tt 2 2? w m * +k l tf+ k i (1 - x)x n ' [ } 



J 



■a M 2 



l-x 



f dxdkj (-l) J+f K + fej)/M 2 

7 16tt 2 ^ M2 _ !!i+a _ £t*l (l-^ 2 ' l J 

X 1 — X 

This form matches that of the equivalent eigenvalue problem in Yukawa the- 
ory [13], with M used as the mass scale instead of /xo and with the replacements 
g 2 — > 2e 2 and I\ — > —1I\. The integrals Jo and J are equal. The solution to 
the eigenvalue problem can be transcribed from [13]; it is 

2M(m 1 -m )(2/iT/ )' ^ M^m/ ' l J 



with z determined by normalization. 

We now must discuss the problem to which we alluded earlier: the appearance 
of new singularities. Since itlq < 1, and for the parameter values of interest 
m > — 1, there will be a pole in the integrand for the % — I — term in the 
equation. The effect looks like a threshold, but there is no state into which 
the system can decay. The fact that the pole can exist is due to the indefinite 
metric. We believe that it is an artifact, and that the correct procedure is to 
define the singularity using the principal-value prescription. That is what we 
shall do for the eigenvalue equation, but the same singularity occurs in the 
normalization integral and in the equation for the magnetic moment. In those 
cases it is a double pole, and we cannot define it as a principal value. We 
shall return to this point presently. Interpreting the singularity in (61) as a 
principal value, we can perform the integrations and then take the limits of 
the PV mass mi — > oo and the physical photon mass /io — > 0. The complete 
result is quite complicated, but for the parameter values of interest we can 



19 



find several much simpler approximate expressions for mo. One is given by 

a l + 21n[// 2 ] 



m Rj 1 — 6 



47r2-|f + gln[/x? 



(65) 



Here we set the physical electron mass M to 1. For all except the largest values 
of Hi in the range of interest, we can just use 



m^l-^(l + 21n[^]). 



(66) 



We must now calculate the integral that appears in the normalization condi- 
tion; it is given by 



a 
2^ 



dx dk 2 ± £(-l) i+J '(l 



X 



hj 



(m 2 — 4m mjX + rn^x 2 ) + k\ 
(x(l — x) — m 2 {l — x) — jifx — k±) 2 



.(67) 



The % — j — term of this integral contains the double pole mentioned 
earlier. Using the prescription discussed in the Introduction, we will define 
the normalization integral as 



a 

2^ 

x 



lim — / dx dz I (1 — x) (ml(l — Ax + x 2 ) + z 



[ml(l — x) + HqX — x(l — x) + z — e] 



mo(l — x) + fi 2 x — x(l — x) + 2 + e]- 

(m^ — AmorrijX + m^x 2 ) + z 



+ £ (-i)^(i -x) 



(x(l — x) — m 2 (l — x) — /i|x — z) 2 



(6J 



with 2; = k 2 . 



This is the closest we can come to defining the double pole as the derivative of 
a simple pole. The major requirement of the prescription for the double pole 
is that, combined with the principal-value prescription for the single pole, 
it should respect the Ward identity. When, we calculate the magnetic mo- 
ment below, we shall also find a double pole in an integral and shall specify a 
meaning for it in a similar way. The double-pole prescription does affect the 
calculation of the magnetic moment; nevertheless, the prescription allows the 
calculation to proceed along lines closely parallel to those of perturbation the- 
ory, and the value of the anomalous moment is consistent with perturbation 
theory. We therefore believe that the prescriptions we have given are consis- 
tent, but we have not explicitly verified that they respect the Ward identity. 
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An interesting consequence of the double-pole prescription is that the normal- 
ization integral is now infrared finite: after performing the integrations we can 
take the limit /xo — > 0. We did not expect this, since in perturbation theory 
the wave function renormalization constant is infrared divergent. The differ- 
ence can be traced to the fact that, with the double-pole prescription, the true 
singularity in the normalization integral is at M = m + \x. Since perturbation 
theory is an expansion around M — m, the singularity in perturbation theory 
is always at \x = 0. In the nonperturbative calculation, the integral involving 
the physical fields has /xo = ^ M — mo, while the integrals involving the PV 
fields do not encounter the double pole at all. The ability to take the mass 
of the physical photon to zero is a useful advantage in performing an approx- 
imate calculation of the electron's magnetic moment because the anomalous 
moment is quite sensitive to a nonzero photon mass [7]. 

Using the double-pole prescription, the normalization condition is quite com- 
plicated. An approximation sufficient for our needs is 



1 + f" ( I (l + 12r^-4mg + 2mg(-7 + 2mg)log(-^_) (69) 




— 2 log(l — m 



We now have the wave function, the bare mass, mo, and the wave function 
normalization as functions of \x\. We can therefore calculate the anomalous 
moment using the wave function overlap formula of [20]. We find 



a [ dx£k m ° 



ir 2 N 2 J x 



1 



m 0+ fc ± _ _^j\2 (I _ m 0+ fc I _ M?+fcj \ 2 

X 1—x' V x 1—x ' 



.(70) 



Here again we encounter a double pole in an integrand, and we use the same 
prescription as above. With that prescription, and making the approximation 
in the second term that itlq = 1, which is good enough for our needs, we find 
that in units of the Schwinger term, £-, the anomalous moment is 



-1 + 2m 2 + 2m 2 {l- m 2 ) ln[- °—} ( 71 ) 

1 TJln 





_ 8 - 6/X! 2 - 3/xx 4 + 2/X! 6 - 6 (-1 + /x x 2 ) 2 log(-l + x*i 2 ) x 1 
( 3(-2 + /X! 2 ) 4 J 'iV 2 ' 

For values of \x\ larger than about 10 this is well approximated by 

m 2 2x1 



a e = m 



(-l + 2^ + 2 m §(l-^)ln [r -^l - )^. (72) 
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With this expression we find that the anomalous moment a e is 1.02 at \x\ = 3, 
1.09 at Hi = 10, 1.13 at ^ = 100 and 1.14 at Hi = 1000. We show a plot of 
the function in Fig. 2. 
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Fig. 2. The anomalous moment of the electron in units of the Schwinger term {-^-) 
plotted versus the PV photon mass, //i. 

Should we view these results as satisfactory or not and what should we choose 
for fj,i7 In the subspace which we have kept, we represent all the processes that 
contribute to the Schwinger term. We also include processes which contribute 
to higher order including all orders. If all the contributions were positive we 
might therefore expect to do better than the Schwinger term, but that is not 
the case: there is a considerable cancellation among the terms which contribute 
to the next order contribution, the Sommerfield-Petermann order a 2 term [24]. 
Thus the best we can expect at this level of truncation is to get an answer 
close to the Schwinger term. If we choose a value of fj,\ anywhere between 
3 times the electron mass and 1000 times the electron mass, we obtain an 
answer within 15% of the Schwinger term. We consider that to be entirely 
satisfactory. 

In order to estimate an optimal value for /ii, we apply the procedure suggested 
in the Introduction: we estimate the error associated with having /xi too small 

2 

as ~ ^f-. We estimate the error associated with having the PV mass too large 

Mi to to 

by performing a first-order perturbation calculation using the projection of P~ 
onto all of the excluded sectors of the representation space as the perturbing 
operator. The calculation proved to be quite challenging. The numerical effort 
to get even an approximate value for the magnitude of the projection of the 
perturbed wave function onto the three-body sector was greater than that 
required for the nonperturbative calculation of the magnetic moment. We 
have done the calculation for two values of fi\, and we find that 



(*'+!*+>! 



+ / |/ii=3 



2 x 10" 



a y 
4W 



(73) 



22 



(<I>' + |^)U =1 oo^8xl0 5 (^J (71) 



We can interpolate between these two values using either linear or logarithmic 
interpolation (asymptotically it should be logarithmic, but we do not know if 
we are in that region); for the present case there is little difference since the 
result is so near Hi = 3. Setting the two types of error equal to each other, 
we estimate the optimum value for \x\ to be between 3.5 and 4. This gives 
us an estimate of the electron's magnetic moment of about 1.02. There are 
unknown factors on both sides of the relation we used to estimate the optimum 
value of Hi, so a considerable uncertainty must be attached to our result. It 
is interesting that the estimate suggests that we should use a value near the 
lower end of the range we show in Fig. 2. In the present case that is the region 
where our result is the best; but we do not know if this is fortuitous. We feel 
that the main points are that the estimate is reasonable and the final result 
is not particularly sensitive to the choice we make, even if we change it by an 
order of magnitude or more. 



4 Light-Cone Gauge 



In this section the notation that we use for light-cone coordinates is 

x ± = -^(x°±x 3 ), x ± = (x\x 2 ). (75) 

v2 



The time coordinate is x + , and the dot product of two four- vectors is 

p ■ x = (p + x~ + p~x + ) — p± • x±. (76) 

The 7 matrices are chosen as follows: 





7* = i I'"" ° I ■ (77) 

a k 



Now we shall again attempt the calculation of the magnetic moment in light- 
cone gauge, but this time by applying the method of regularization and renor- 
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malization that gives a light-cone perturbation theory equivalent to the stan- 
dard, dimensionally regularized Feynman series to all orders in the coupling 
constant [16]. This method includes the introduction of higher derivatives into 
the Lagrangian, along with the use of three PV electrons and one "technical" 
photon. Certain cutoffs, discussed below, are also required. From the perturba- 
tive equivalence with Feynman methods, it is clear that the method preserves 
gauge invariance in the limit where regularization is removed, at least at the 
level of perturbation theory. The Fock-space truncation will break gauge in- 
variance, as it does in the Feynman gauge. The differences between the results 
in this section and those in the previous section give some measure of the 
effect of the breaking of gauge invariance by the truncation. 

The starting form of the Lagrangian density is 



'=-!e(-^4 + m 




4^ '' ""M A* A* 

+ E ~# Wd, - rm) $ + e^YA^', (78) 

1=0 u i 

where F mj = d^A v - d v A^ A^ = A^ + A^, A_ d = 0, V' = £?=o$> 
E?=o^ = E?=o^™i = £f=o*W = 0, i/o = 1, and d± = (d ± d^/y/2. 
We should comment on the photon fields labelled with j — 1, which we call 
the technical fields. In [16], it was explained that these fields implement the 
Mandelstam-Leibbrandt (ML) prescription for the spurious singularity. In the 
nonabelian theory, it is known that this prescription is necessary for a consis- 
tent perturbation theory. In the abelian theory, it is not known whether the 
ML prescription is necessary or if, for instance, the principal-value prescription 
will suffice [25] . On that basis, one might wonder whether the j '• = 1 fields are 
needed in QED. However, the j '• = 1 fields are necessary for the proof in [16] 
of the equivalence between light-cone and covariant Feynman formulations of 
perturbation theory in either the abelian or the nonabelian case. What we 
find in the present calculations is that the j '• = 1 fields are necessary even for 
this relatively simple problem in QED. 

For the UV regularization of the photon field, we use higher derivatives, which 
do not break gauge invariance in QED. To construct a Hamiltonian formal- 
ism on the light-cone, we rewrite the system in a form which includes only 
first derivatives of the fields with respect to x + . To regulate the theory in a 
way which guarantees the equivalence between light-cone and Feynman per- 
turbation theory, we must also introduce PV electrons with flavor-changing 
currents. This breaks gauge invariance, as it does in the Feynman gauge. How- 
ever, as in Feynman gauge, that breaking of gauge invariance does not require 
counterterms, since we can take the limit m li2 ,3 — > oo. The truncation of 
the Fock space excludes fermion pair creation, so there are no fermion loops; 
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therefore, after taking the limit e — > (i.e. removing the light-cone regular- 
ization \k-\ > s > 0), we can take the limit 777,1,2,3 — > oo at fixed A. We then 
find that gauge invariance is restored without using counterterms in the La- 
grangian. The calculations that we carry out here confirm that these limits 
indeed turn out to be finite. If we had included enough of the representation 
space to allow pair production to occur, then we would not be able to take 
the limit 7711,2,3 ~~ y 00. In that case, we would have to add counterterms and 
also renormalize the electric charge e. 

To get a canonical light-cone formulation and determine the Hamiltonian P~, 
we decompose the bi-spinors ■*/>' into two-spinor components, ip[ = {ipi,+, fy,-}, 
and express the components ■?/>;, _ in terms of other field variables using the 
canonical constraint equation on the light-cone. As shown in [16], it is con- 
venient to use <fj = d^Aj and ip = d^A^ in place of the A + j as degrees of 
freedom, and we shall make that choice. With all these choices, P _ can be 
written as 



p-= [dx- ((fx 1 - W A2( ir 
J J AiU (1 + ^1 



sH -HB 



+ 



fc=i, 

:-iv +:l 



-4,, + 4fl + CV-4 



k ^J?\ L ^^l\ ^3 



E A ^i 1 



'_]V'+i 



fc=l,2 



d 2 ± \ 

J2 ) A *>i + — 



V/(i-^)^ 



+^E^rV^(5i-m2)9i 1 ^+-ev / 2^ + r 1 V? + E 9 k A k 






fc=l,2 



(79) 



where n^j 



<5L 



<5(a+A feiJ ) 



are canonical momenta conjugate to Akj, and we have 



defined A± = Ai<j 1 + A 2 <J2 and d± = (c\, <9 2 ). 

We decompose our field variables in terms of creation and annihilation oper- 
ators acting in light-cone Fock space: 



A kd (x) = (2tt)- 3 / 2 d 3 P x: 



a jC,k{p) e %PX + ^-C 



JO) 



^ , hr>+ U 



2P + ( ^ " 1 



n fcj (a:) = -2(27r)- 3 / 2 /rf 3 p5: * 

C=o,i > 



^ ^4 - lV (a j<k (p)e- ipx - /i.e.) , (81) 



2 V^ 2 
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^(x)=i(2n)- 3 / 2 A 1 + 



A 



2\ -J'/2 



rf 3 p 



a j(k(p) e ipx — h.c. 



where s is a spin projection. The integrals are defined as 

oo 

d 3 p = / d 2 p ± / dp + , 



82) 



i/>i,+,s(x) = 2- 1 / 4 (2tt)- 3 / 2 />p (6 M (p)e-^ + < s (p)e^) , (83) 



34) 



where t> and £ > are regulation parameters which must be taken to zero 
in a prescribed order, as was shown in [16], and will be discussed below. The 
nonzero commutation relations for the creation and annihilation operators 
have the following form: 



[ a jCk(p), a fC'k'(p')] = (~ 1 ) <:s kk'S jj '5 cc 8 3 (p - p'), 
[ai(p),aJ(f/)] = (-l) i+1 ^« S (p-p'), 

hs(p),K,Ap')] = ldiAp)i4,Ap')} = vi5u>6 ss ,5 3 (p-p') 



35) 



Substituting these decompositions into (79), and keeping only the terms nec- 
essary for our approximate calculation of the anomalous magnetic moment, 
we obtain 

l,s j 

+ Y:Jd 3 k(-l)<E K (k)a+ cx (k)a KX (k) + j^ J2 Jd 3 qjd 3 k 



J'.CA 1 



A 



x (aj(k)K.{Q + k)b v , a (q) + a+{k)b+ s (q - k)k'Mj v2 J/2 



(i+r^ 1 



E 



CA ■ hk+ f4-i 



X 



g(-A) 

g+ 



^ -A/2 + 



(q + fe) ( _ A ) „ /.' 



(g + A;)- 

'9(A) 



^ 



s,A/2 



(-A) 






(A) 



$s,-\/2 



+«; C A(*)6f.(9 - fc )^^9) ^^>A/ 2 + (r/ _ A . )+ 

+ 7f ( (g + fc)+ ~ g^J ajCA(A:) ^ s(g + k ) b i',-s(qK-x/2 



K (A) 

A;+ 
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+- 



/ m t 



m v 



a Kx( k )Ks(l ~ k)b v ^ s (q)5 s ^ /2 



V2\(q-k)+ ^ 
where we use helicity components for the operators ctj^ip)'- 



a iCA = —ft ( a j(i + *Aa iC2 ) , 



A = ±l, 



36) 



J7) 



and we have defined 



P(A) 



—f= (pi + i\p 2 ) , P 2 ± = p\ + p\ = E P(A)P(-A), 
V 2 A=±l 



e M , EL+f*, * jW , -4^. *«<*) - kl + CA ' 2 



2^ 



2k+ 1 



A 2\J 



2A;+ 1 + C 



A 2\J 



Now we truncate the Fock space and write the wave function as 



\p) = (e fi(p)Ks(p) + E / rf3 <? fiM <i)Ks(<i)4(p - ?) 

\ l,s l,j,s 

+ E /A/5cA(P.?W.(?)a; cA (p-?) > )|0>. (89) 

As in the previous section, we can choose a state with p± = 0. The eigenvalue 
problem in this subspace is 



M 2 
P+p-\p) = —\p). 



(90) 



Acting with P , as given by (86), on the state (89), Eq. (90) gives the following 
relations among the coefficients of the basic state vectors. For b^~ s (p)\0) we have 
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where x = q + , z = (p — q)\. For 6 Z + S (g)a^(p — g)|0) there is 



0, (91) 
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For bf {q)aj^ x ip — q)\0) we find 
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Expressing the wave functions, ffj(p,q) and //^(p, <?), in terms of the /f(p) 
and substituting into the relation (91), we obtain the eigenvalue equation 
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(94) 



where we have defined z = q\ . 
We write this equation in the form 



ml-M 2 = -{AM 2 + 2Bm + C) + ^(AC - B 2 ) 
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(95) 
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where A, B, and C are real functions of itiq, M, and the regularization pa- 
rameter A. To restore Lorentz invariance and gauge invariance, except for the 
breaking due to the truncation, we must now take the following limits in the 
order given [16]: 

e — > 0, /j,—*0, v — > 0, mi, m2,m3 — * oo, (96) 



We can now find the value of mo by setting M equal to 1 and solving (95) nu- 
merically. We can find the leading order solution by expanding the expressions 
for A, B, and C for large values of A. We find 

A = ^\nA 2 + ..., B = -^\nA 2 + ..., C = ~- In A 2 + . . . .(97) 



Keeping only these terms, we can solve Eq. (95) and obtain 
1- if- hi A 2 ) _ A(£i n A 2 ) 2 



H-(flnA 2 )+i(fmA 2 )' 



This leads to the following expression in the lowest order in coupling constant: 
m 2 = M 2 (l-|hnA 2 V (99) 



This expression agrees with Eq. (66) (and with the known one-loop result for 
the self-energy of the electron). As expected in QED, m\ < M 2 . Therefore, the 
integrand of equation (94) does contain a simple pole that we have integrated 
using the principal-value prescription, as we proposed in the Introduction. 

For either value of s, we can parameterize the constants, ff, in the wave 
function in terms of A and B as 
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To obtain a normalized wave function we must divide by the norm N, which 
is given by the expression 
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(101) 



Here again we encounter the double pole, which we treat according to the 
prescription (4) (see the Introduction). 

The result of calculating the integrals in (101), and taking the proper limits 
specified in (96), is a very complicated expression that is a function of mo, 
M, and A. The limits must be taken only after performing the integrals of 
Eq. (101). 

To calculate the anomalous magnetic moment of the electron, we need to 
use the norm of the eigenstate (101) and calculate the matrix element of the 
current component j + {x) between the states (89). Using the method of [20], 
we can write the expression for the anomalous magnetic moment a e as 



2M d . . , . .. 

■j^i-Q—{P± = Pi; s\r{0)\p± = 0; -s) | P1=0 , 



(102) 



where the relevant part of the current (containing only b^ s and 6+ s ) is 



J + (0)= E j d 3 qid 3 q 2 bt iyS ( qi )b hyS (q2). 



(103) 
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Using the expression in (89) and Eqs. (92) and (93), we get 
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In the calculation of this integral, we have again used the prescription (4) for 
the double pole. 

One can show that after removing the regularization, (96), we have f( = f£ = 
fi = ffo = 0- Then the calculation gives an expression which can be written 
at large A as 
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(105) 



Substituting the values of mo, f , and N, found from (95), (100), and (101) at 
fixed a and A, and setting M — 1, we get the value of a e . At lowest order in 
a (when m — > M, N — > / ), the expression (105) gives the Schwinger term 



Lvp — 

c 2-K 



In Fig. 3 we show the value of the anomalous magnetic moment, a e , in units 
of the Schwinger term (^), plotted against the UV regularization parameter 
A. 




1000 



Fig. 3. The result of a light-cone calculation of the anomalous moment of the 
electron in units of the Schwinger term (tP-) plotted versus the UV regularization 
parameter A. 
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5 Discussion 



In this paper we have presented nonperturbative calculations of the electron's 
magnetic moment. In each case we used the light-cone representation and 
solved the eigenvalue equation in a basis truncated to include only states with 
one electron and states with one electron and one photon. Once we found 
the one-electron eigenstate in the truncated subspace, we then calculated the 
magnetic moment from the eigenstate. 

We found that there were three problems that we had to overcome in order 
to make a successful calculation: the problem of uncancelled divergences, the 
problem of new divergences, and the problem of maintaining gauge invariance. 
Our solution to the problem of uncancelled divergences is to keep the PV 
photon mass (or the, nearly equivalent, UV regulator A in the light-cone gauge 
case) finite. For that method to succeed we had to find a large range of PV 
photon masses for which the error due to including the negative metric states 
in our wave function and the error due to the truncation of the representation 
space were both small. We did find such a range. If one accepts the errors 
shown in Figs. 2 and 3 as being not so large as to render the calculation 
useless, we must only estimate the optimum PV mass as lying between three 
electron masses and one thousand electron masses in order to perform a useful 
calculation. The estimate we obtained, ~ 4m e , not only lies in this range but, 
perhaps fortuitously, lies in the part of the range in which our answer is nearest 
to the correct answer. 

In our calculations we encounter integrals whose integrands contain poles, 
which we define as principal values, and double poles. In the case of the dou- 
ble pole we have provided a prescription which is the nearest we can come 
to defining the double pole as the derivative of a principal value. With these 
prescriptions we obtain a successful calculation. The ultimate test of the pre- 
scriptions is whether or not they respect the Ward identity. We have not 
proven that the prescriptions preserve the Ward identity beyond the present 
calculations. An unexpected feature of the prescriptions is that the wave func- 
tion normalization is now infrared finite, whereas it is infrared divergent in 
perturbation theory. We have therefore been able to carry out the present 
calculations using zero for the mass of the physical photon. 

The most difficult problem we had to solve was the problem of maintaining 
gauge invariance. If the standard light-cone procedures are applied, regulated 
with any combination of PV fields that we tried (or, regulated with a mo- 
mentum cutoff), the result is not a successful calculation. The failure of the 
standard calculation to produce useful results can be traced to the loss of 
gauge invariance. Understanding the details of that loss of gauge invariance is 
an interesting unsolved problem. We have produced successful calculations by 
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the use of Feynman gauge and by the use of light-cone gauge regulated with 
both higher derivatives and PV fields. It should be noted that in either of 
the formulations which resulted in a successful calculation, A~ is a degree of 
freedom, whereas in the unsuccessful light-cone gauge calculation A~ satisfies 
a differential constraint equation. 

An important observation is that, with the use of the PV fields, the oper- 
ator P~ can be constructed without inverting any covariant derivatives. In 
the past, most calculations for gauge theories in the light-cone representation 
have been done in light-cone gauge due to the apparent need to invert a co- 
variant derivative in any other gauge. With the use of appropriately coupled 
PV fields, that impediment to the use of other gauges, including covariant 
gauges, is removed. In the present work the calculations in Feynman gauge 
are considerably simpler than the successful calculations in light-cone gauge. 
Looking forward to the nonabelian case, it is not clear that this will remain the 
case. Use of the Feynman gauge would involve the inclusion of Faddeev-Popov 
fields, whereas presumably that would not be necessary in light-cone gauge. 
Also, the nonabelian formulation corresponding to Sec. 4 has been shown to 
give perturbative equivalence to covariant methods to all orders in perturba- 
tion theory. No such demonstration has yet been given for other gauges. 

Even when we have a formulation which, in the absence of truncation, would 
preserve gauge invariance, the truncation will break gauge invariance. We have 
argued that the question of whether or not that breaking of gauge invariance 
is acceptable is more a question of accuracy than symmetry: if our broken 
answer is close to the correct, unbroken answer without truncation, it will be 
a useful answer. The difference between the answers we obtained here in the 
Feynman gauge and in the light-cone gauge provides some test of that idea. 
The answers differ by some two percent. We believe that the fact that they 
are close to each other is due to the fact that they are close to the answer 
which would result from a calculation without truncation. 

The objective of the calculations presented here was not to improve on the 
very accurate calculations of perturbation theory. Rather, the objective was 
to verify that the nonperturbative methods we have developed really do repre- 
sent a nonperturbative approximation to QED. We have succeeded in this, not 
only because the answers we get are as close as can be expected to the correct 
answer but also because an expansion of our approximate, nonperturbative 
answer in powers of the coupling constant reproduces the finite terms of stan- 
dard QED perturbation theory, through the order allowed by our truncation 
of the representation space. We get additional terms up to all orders in the 
coupling constant, which can be identified as some of the terms of standard 
perturbation theory. Therefore, it is reasonable to hope that our nonperturba- 
tive methods will produce useful calculations for problems which perturbation 
theory cannot address, such as hadron bound states. 
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